Cavity-induced temperature control of a two-level system 
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We consider a two-level atom interacting with a single mode of the electromagnetic field in a 
cavity within the Jaynes-Cummings model. Initially, the atom is thermal while the cavity is in a 
coherent state. The atom interacts with the cavity field for a fixed time. After removing the atom 
from the cavity and applying a laser pulse the atom will be in a thermal state again. Depending on 
the interaction time with the cavity field the final temperature can be varied over a large range. We 
discuss how this method can be used to cool the internal degrees of freedom of atoms and create 
heat baths suitable for studying thermodynamics at the nanoscale. 

PACS numbers: 42.50.Pq, 32.80.Pj, 05.30.-d 



The Jaynes-Cummings model [lj (JCM) is a simple but 
powerful model describing the interaction between a two- 
level atom and a single mode of the radiation field. While 
being exactly solvable it offers a large variety of genuinely 
quantum phenomena like collapses and revivals in the 



inversion of the atom 
experimentally as well 



which have been observed 



For a field prepared in a coherent state, the state of the 
atom will be almost pure at half of the revival time if the 
atom is initially in a pure state @, 0|- However, a more 
realistic model would involve a thermal initial state for 
the atom. The thermal contribution to the initial state 
of the field may be neglected as long as the number of 
coherent photons is sufficiently larger than the number 
of thermal photons 0]. Using thermal states allows for 
an investigation of the thermal properties of the JCM, 
i.e., its applicability for problems like the initial state 
preparation in quantum computing 0] , cooling of atoms 
[lfj , [ill H3 |. or implementation of quantum thermody- 
namic machines 



impli 



In the following we will discuss a model where an atom 
in a thermal state enters a cavity prepared in a coher- 
ent state. By obtaining a closed form for the reduced 
density matrix for the atom we will show that after the 
collapse the state of the atom is independent of its initial 
state. After a fixed interaction time the atom is taken to 
leave the cavity and to interact with a laser field, which 
is treated as a semi-classical driver. For an appropriate 
laser field the final state will be thermal as well. De- 
pending on the interaction time with the cavity the final 
temperature can be varied over a large range, leading 
to cooling or heating of the atom. We will present an 
expression for the minimum and maximum temperature 
that can be achieved. Finally, we will discuss applica- 
tions of the method to cooling of the internal degrees of 
freedom of atoms and creating heat baths suitable for 
studying thermodynamics at the nanoscale. The whole 
procedure of our proposal is depicted in Fig. [TJ 




Laser 



FIG. 1: Outline of the procedure: A two-level atom in a 
thermal state (occupation probabilities represented by black 
dots) interacts with a cavity prepared in a coherent state \a). 
After a time t the atom leaves the cavity. A laser pulse is 
applied to the system, resulting in a thermal state with a 
different temperature. 



The total system is described by the Hamiltonian 

H = Ha + Hf + Hj, (1) 
where the atomic Hamiltonian Ha is given by 

(2) 



with AE being the energy splitting. The field Hamilto- 
nian Hp is 



Hp = huj 1 a a + — 



it; 



(3) 



with u> being the frequency of the single mode and a 
being the annihilation operator of the field. Being in 
resonance, we have w = AE/h. Using the electric dipole 
and rotating wave approximation the JCM interaction 
Hamiltonian is given by 



Hi 



g& + a + g*a a\ 



(4) 



where g is the coupling constant for the atom-field inter- 
action and & are the atomic transition operators. Par- 
ticularly, the coupling constant is given by 



g = d 



he V 



(5) 
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where d is the atomic electric dipole matrix element and 
V is the mode volume. 
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We restrict ourselves to the field being initially in a 
coherent state |a) and the atom being in a thermal state 
described by the density operator 

p A (0) = zT 1 exp{-pH A ) = Pe(0)\e) (e| + [1 -p e {0)]\g) (g\, 

(6) 

with Z being the partition function, fi the inverse tem- 
perature, p e the probability to find the atom in its excited 
state |e), and \g) its ground state. 

The time evolution of the full system is then given by 

p(t) = Pe (0)U\e,a)(e,a\fr + [l-p e (Q)]U\g,a)(g,a\& 

= Po (0) (*)) (V-e (t) | + [1 - Pe (0)] | % (t)) <V> 9 (t) | , 

(7) 

where [/ is the time evolution operator of the full system. 

In order to obtain the effective time evolution for the 
atom alone, the degrees of freedom corresponding to the 
field have to be traced out jl4j |. Here, the partial trace 
over the field is given by 

p A {t) = Tr F {m))(m\} = 5>M*)>W*)|n>. (8) 

n 

Since pA is Hermitian and has unit trace, the atom is 
effectively described by the diagonal element pn and the 
off-diagonal element poi- 

In the following we first consider the case where the 
initial state is |e,a). Then, the full time evolution is 
given by (see, e.g., [II]) 

\M*)) = E (e- in " +lt/2 |+« + i) - e*W/ 2 K+1 ; 



V2 Vri. c 



-i(ra+l/2)wt 



(9) 



where the n-photon Rabi frequency fl n = gyfn and the 
n-photon eigenstates of the atom-field system, 

\± n ) = -j=(\g,n+l)±\e,n)), (10) 

have been used. The time evolution of the reduced den- 
sity matrix element p\\ before the revival time have been 
studied extensively (see, e.g., [Ill) and is given by 



, s 1 1 s ft 2 

Putt) = - + - cos(2 5 t) exp y-^T i 



(11) 



with tc being the collapse time tq = v2/ 1 g. For an atom 
initially in \g) the result is 



/ N 1 1 ic\ \ ft 

pu(t) = -- -cos(2 3 <)cxp y-^J 



(12) 



Therefore, after the collapse the diagonal elements are 
constant and pa = 1/2. 

The off-diagonal element poi (again, first for the atom 
initially in |e)) is given by 

Pox(t) =^2(ip e (t)\n,g)(n,e\i; e {t)). (13) 



Evaluating the summands p^ using Eq. (J9j> and Eq. (fTO]) 
leads to 



Poi (*) = 



-\LUt 



2a* 



(O n _|_i + Q n ) 



(^n+l — 



(14) 



with w(n) being the Poisson distribution. The first term 
inside the square brackets oscillates at a much higher fre- 
quency than the second and results only in a random 
phase, which vanishes after summation. In the high- 
photon limit y/n may be approximated by (see @) 



(15) 



Analogously, the difference of the Rabi frequencies can 
be expressed as 



fi 



n+l 



fir, 



2g(Vn + l 



1 



2y/E 8 



(16) 
(17) 



Plugging only the leading order into Eq. (TT4"1) and replac- 
ing the sum in Eq. (|13[) by an integral over a Gaussian 
distribution leads to 



poi(t) = --exp[i(o;i 



■ (j))] sin 



gt 

2 s/n 



(18) 



where is the initial phase of the radiation field. Us- 
ing the same approximations for the atom initially in its 
ground state yields the same result for pm(t). Therefore, 
after the collapse the atom evolves totally independent 
from its initial state. A comparison of Eq. (fTS]) with the 
numerical solution of the full time-dependent Schrodingcr 
equation is shown in Fig. [2] Apart from the collapse and 
revival phase there is excellent agreement. This further 
shows that the random phase approximation applied to 
Eq. (fTl)) was perfectly justified. 

Since the diagonal elements of pa are both at i the 
Bloch vector only moves within x — y plane of the Bloch 
sphere. Furthermore, the phase oscillates at to, i.e., in 
the rotating frame only the initial phase <fi is relevant [see 
Eq. (|18p]. Therefore, in order to obtain a thermal state 
one always has to apply a n/2 pulse to the system (see 
Fig. El) ■ Since the pulse diagonalizes pa, the probability 
to find the atom in its excited state after the pulse p e (t) 
is given by the smallest eigenvalue of pa- Computation 
of p e (t) yields 

Pe(t)4(l-sin^=). (19) 
This can also be expressed as a temperature using 

T= (20) 
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FIG. 2: Comparison of the real and imaginary part of Eq. (|18[) 
(solid lines) and the solution of the full time-dependent 
Schrodinger equation. Initial states for the atom were \g) 
(Repoi: crosses, Impoi: boxes) and |e) (Repoi: diamonds, 
Impoi: triangles), (n = 36, g = AE, and = 0) 




FIG. 3: Illustration of the 7r/2 pulse acting on the Bloch 
vector of the atom. 



This temperature should be considered as a param- 
eter characterizing the mixedness of the output state 
rather then as an indication for a stable thermal state 
proper. Quantum objects prepared like this (or ensem- 
bles thereof) would then constitute resources for further 
applications like quantum information processing. 

Equation (p~9|) suggests that at half of the revival time 
the atom will be in its ground state (i.e., T = 0). How- 
ever, this minimum temperature would only be reached 
for infinitely large n, for which it would take an infinitely 
long time to reach this state. In order to determine the 
actual minimum temperature a correction for finite n is 
required. A correction to Eq. (|19p can be obtained by 
including the next order in Eq. (fTT)) . Close to half of 
the revival time the sine in Eq. (|14p is near its maxi- 
mum and can be approximated by a second order Taylor 
expansion, which leads to a final result of 



y V 2 / 32n 



(21) 



Using the next order in Eq. (|15p as well leads to an 
additional correction in 0(l/n 2 ). Putting this p e into Eq. 
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FIG. 4: Minimum temperature T m i n over average photon 
number n. 



pop gives the minimum temperature T m j n (fi) as shown 
in Fig. |H Temperatures as low as 0.2 AE/ks can be 
obtained, which correspond to an occupation probability 
of the excited state of the order of 10~ 3 . 

In order to determine the maximum temperature that 
can be reached we require that the collapse must have 
taken place [i.c,. the difference in the occupation proba- 
bilities Eqs. (HI|) and (12]) is negligible compared to the 
difference induced by the laser] . Requiring the former to 
be smaller by a factor of 10, this can be expressed as 



, t 2 \ gt 
10 cos(2g/j) exp [ j J = sm 



'C 



2 VI' 



(22) 



The cosine on the left hand side may be replaced by unity 
without violating the above requirement. For large Ti- 
the right hand side can be approximated linearly in t, 
resulting in 



t 2 

10 exp [ — j 



gt 



(23) 



Solving for the appropriate cavity interaction time t and 
using Eqs. (fl"9|) and (|20|) leads to a maximum tempera- 



ture T max of 



AE 



4:y/n-\- n) 



(24) 



k B log 



y/W(400n) 



where W(-) denotes the Lambert W function, i.e., the 
inverse function of f(x) — xe x . Figure [5] shows the de- 
pendence of T max on n. 

These results show that the temperature can be tuned 
over a large range, which depends only on the average 
photon number n, the coupling time t, and the energy 
splitting AE. However, there are some other applications 
for this procedure, which are realizable within present ex- 
perimental setups. A rather obvious one is the cooling 
of the internal degrees of freedom of atoms. However, an 
implementation using a cavity would be extremely diffi- 
cult as the frequencies relevant for cooling are in the MHz 
range, where the coupling constant g is much too small to 
observe any effects [due to the u> dependence in Eq. J5J)]. 
A much more promising implementation could be realized 
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FIG. 5: Maximum temperature T max over average photon 
number n. 

using circuit quantum electrodynamics (QED) [l6l.[l7l]. in 
which the atom is replaced by a Cooper-pair box and the 
cavity is implemented by a one-dimensional resonator. 
There, the coupling constant is sufficiently large even in 
the relevant frequency range. Although dephasing plays 
an important role in circuit QED, the dephasing time is of 
the order of several hundred nanoseconds [18| |. while the 
required time for reaching the minimum temperature is 
about one order of magnitude smaller for experimentally 
feasible parameters. Using our procedure thus might 
lead to lower temperatures than currently employed tech- 
niques. Besides circuit QED, other implementations in- 
volving a Jaynes-Cummings Hamiltonian with a tunable 
coupling constant may prove useful as well. 

Another interesting application of this procedure could 



be the realization of tiny local baths. Local baths are an 
important ingredient in nonequilibrium quantum ther- 
modynamics [13| , where it is necessary to create and con- 
trol a temperature gradient on a nanoscopic scale. This 
could be used to investigate transport behavior [2(| 
or quantum thermodynamic machines [2l[. Using our 
framework to repeatedly set a temperature of a single 
two-level system could act as such a local bath as long 
as the cavity is reset after each step and the tempera- 
ture control happens on a much smaller timescale than 
the other processes within the system (i.e., strong bath 
coupling) . 

In summary, we have shown that the temperature of 
a two-level atom could be efficiently controlled via a res- 
onant interaction with a cavity. Depending only on the 
interaction time with the cavity, it should be possible to 
tune the final temperature over a large range. The ex- 
pression for the reduced density matrix of the atom has 
been obtained in the high photon limit using a system- 
atic series expansion and has been verified by comparison 
with the solution of the full time-dependent Schrodinger 
equation. Besides temperature control our procedure 
may prove useful for cooling various microscopic systems 
or realizing local baths in nanothermodynamics. 

We thank T. Pfau, M. Michel, M. Henrich, F. Rempp, 
G. Reuther, H. Schmidt, H. Schroder, and M. Youssef for 
fruitful discussions. 
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